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ABSTRACT 

Cosmological parameters affect observables in physically distinct ways. For example, the baryon density, io b , 
affects the ionization history and also the pressure of the pre-recombination fluid. To investigate the relative 
importance of different physical effects to the determination of LOb, and to test the cosmological model, we 
artificially split u>b into two 'metaparameters': io be which controls the ionization history and u> bp which plays 
the role of uj b for everything else. In our demonstration of the technique we find ui b = .0229 ± .0012 (with no 
parameter splitting), w bp = .0238 ± .0021, uj be = .0150 ± .0034 and u bp -u be = .0088 ± .0039. 
Subject headings: cosmology: theory - cosmology: observation 
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1. INTRODUCTION 



As predicted (So ergell 119951 iKnoxl [l995t Uungman et alJ 
119961) . observations of the cosmic microwa v e background 
(CMB) anisotropies (e.g . iKuoetalJ ( 120041) : iBennettetaD 
(2003); Readhead et al. (2004J)) have provid ed very tight 
constr a ints on cosmological parameters (e.g . ISnergel et al.l 
( 120031) : iGoldstein et alJ (120031) : iRebolo et all J2004l) L These 
constraints are possible because the statistical properties are 
sufficiently rich and, given a mod el, can be ca lculated with 
very high accuracy (e.g. Hu & Dodelson (2002)). 

One must bear in mind though that these determinations are 
highly indirect and model-dependent. It is therefore useful to 
have tools for testing the model, and for gaining better under- 
standing of the particular physical processes important for a 
given constraint. Toward these ends, we explore use of cos- 
mological 'metaparameters'. 

A given parameter, x, may lead to observational conse- 
quences through more than one distinct physical effect. Such 
a parameter can be split into more than one metaparameter, x a , 
x b , x c , ... each of which controls a different physical effect. 
This approach to data analysis has been developed indepen- 
dently and applied recently by Zhang et al. (2003) who call 
it 'parameter splitting.' Their motivation was to marginalize 
over physical effects that could not be calculated with suffi- 
cient accuracy. One can also use the split into metaparame- 
ters to test the model (by checking if x a = x b to within errors), 
and to understand where the constraints are coming from (by 
comparing a(x a ) to a(x b )). 

In this paper we explore one example of a split into cos- 
mological metaparameters. In particular, we split the baryon 
density into two parameters: one that controls the ioniza- 
tion history, u> be , and one that controls the inertia of the pre- 
recombination baryon-photon fluid, uj bp . In section 2 we de- 
scribe the dependence of the angular power spectrum on these 
two variables. In section 3 we briefly describe our calcula- 
tions. In section 4 we show and discuss our results. Finally, 
in section 5 we conclude. 

2. DEPENDENCE OF C L ON u BE AND u BP 

The baryon density affects the evolution of CMB temper- 
ature anisotropies in two distinct ways. First, in the pre- 



recombination plasma, the higher the baryon density, the 
higher the inertia of the fluid due to the baryon's mass. Sec- 
ond, the ratio of baryons to photons determines the history of 
the number density of free electrons (prior to reionization). 

These two different physical effects, one related to the 
baryon's mass and one related to the large Thomson cross sec- 
tion of the electrons that charge-balance the baryons, lead to 
different observational consequences. They therefore lead to 
two different ways to determine the baryon density. 

We can study these effects separately by splitting to b = ft b h 2 
into what we will call ui be and ui bp where the e indicates elec- 
tron and the p stands for pressure. Operationally, we calculate 
the mean number density of electrons, n e , assuming io b = uibe , 
using the recombination routine RECfast ( Seager et al. 1999). 
Then we calculate Q using CMBfast (Seliak & Zaldarriaga 
1996) with uji, = LObp and the n e (z) from RECfast. 

2.1. Dependence on ujbp 

The response of C/ to the two baryon densities can be seen 
in Fig. [2 First we concentrate on the multiple effects of vary- 
ing Lo b p. Decreasing ui bp leads to three different physical ef- 
fects: sound speed increase, acoustic oscillation zero-point 
shift, and reduction in 'baryon drag'. We briefly review the 
effects here though they are all discussed at length in the re- 
view by Hu& Dodelson (200l. 

Adding non-relativistic baryons decreases the sound speed 
of the fluid, given by 
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(1) 



where R = 3p b /p 1 . Reducing uib p decreases R and therefore 
increases the sound speed, causing the oscillation pattern to 
shift slightly to lower /. 

In the absence of baryons, if the gravitational potential "J 
is constant then the effective temperature, 9o + ^, oscillates 
about zero; that is, the competing effects of pressure and grav- 
ity cancel when 6o = —9, Adding baryons reduces the pres- 
sure, meaning the baryons must collapse further into a poten- 
tial well before gravity and pressure balance. The zero point 
shifts to 6 =-(l+R)^> and therefore 6 + * oscillates about 
—R^f. The offset enhances odd peaks (for which the effective 
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FIG. 1. — Dependence of C; on the metaparameters. Solid curve is C\ for 
the fiducial model with u}j,„ = u)}, e = 0.023. Dashed (dotted) curve is Q for 
fiducial model parameters except for w^p = 0.018 (uif, e = 0.018). The fiducial 
model is a typical sample from the MCMC chain of our six-parameter model 
space constrained by the WMAPext dataset. 



temperature is positive in potential hills and therefore is 
a boost) and suppresses even peaks (for which the effective 
temperature is negative in potential wells). Decreasing u>bp 
therefore suppresses odd peaks and enhances even peaks. 

In the above paragraph we assumed a constant which is 
a good approximation for the matter-dominated era, but not 
for radiation domination. For modes that enter during radi- 
ation domination, pressure resists transport of material into 
the potential well, causing the potential to decay as expansion 
dilutes the over-density. Thus for modes that entered early 
during radiation domination, by the time of last-scattering the 
gravitational potential is insignificant and the oscillation is 
about Go = 0; there is no offset to the oscillations and there- 
fore no modulation of the even and odd peak heights. 

At smaller scales, the dominant effect of a reduction of uib p 
is a reduction in 'baryon drag' . The baryon drag effect can be 
thought of as due to the increasing value of R over time. As R 
increases in time, the sound speed decreases, so the oscillation 
frequency decreases. Since energy/frequency of an oscillator 
is an adiabatic invariant, this decrease in frequency is matched 
by a decrease in oscillation amplitude. Since R oc R, decreas- 
ing ojbp decreases the amount of baryon drag, and the power is 
enhanced. Note that the baryon drag effect is is distinct from 
the photon diffusion we discuss in the next subsection. This 
distinction is clear from the fact that there is baryon drag even 
in the tight coupling limit. 

The scale of matter radiation equality projected to today 
comes out at / eq ~ 170 jRnox et alJl2001l) and therefore the 
first peak entered slightly before matter-radiation equality. 
Both effects (zero-point offset and baryon drag) are important 
for the first peak and partially cancel each other out. For the 
second peak, the effects add so that decreasing u>bp enhances 
power. For the third peak, having entered earlier during radi- 
ation domination, the offset effect is sub-dominant to baryon 
drag so power is enhanced. 



2.2. Dependence on u)be 

Even prior to recombination, the mean free path due to 
Thomson scattering is non-zero. Photons thus diffuse, damp- 
ing the power spectrum on small scales. The comoving damp- 
ing length is given by the distance the photons can random 
walk by the time of last scattering. Decreasing u>be decreases 
the number density of baryons, and therefore, for fixed 
ionization fraction, the number of free electrons, 
The resulting increase in mean free path increases the damp- 
ing length. 

The effect of varying uib e on the damping length is compli- 
cated by the fact that decreasing u>b e increases the number of 
photons per baryon and therefore alters x e (z). More photons 
per baryon delays recombination. This means that at a given 
time more atoms are ionized (tending to decrease the damping 
scale), and it also means a delay in recombination giving more 
time for th e random walking (tending to increase the damping 
scale). IHu & White! Jl997h found that the net result is that the 

-1/4 

damping length scale is roughly proportional to uj h . Reduc- 
ing LOj,e therefore increases the damping length, leading to the 
increased suppression for the dotted curve seen in Fig.^ 

The damping length projected from the last-scattering sur- 
face to here, a comoving distance D away, gives rise to an 
angular scale, ir/lp. Below we us e In as a f unc tion of cos- 
mological parameters as written in IHu et alj d2001h with the 
difference that we numerically evaluate the redshift of last- 
scattering, z*, as the peak of the visibility function, and re- 
place LOb With Ulbe- 

3. CALCULATING THE LIKELIHOOD 

What we want to know is, given the data and any other as- 
sumptions we make about the world, what is the probability 
distribution of the parameters? This posterior probability dis- 
tribution can be calculated by use of Bayes' theorem which 
states: 

P(9\d)^P(d\9)P plio[ (9). (2) 

where d refers to data and the proportionality constant is cho- 
sen to ensure J P(9\d)d9 = 1 . With a uniform prior this simply 

reduces to P(9\d) oc P(d\9). This probability of the data given 
the parameters is, when thought of as a function of the param- 
eters, called the likelihood, C(6). 

Often we are interested in the posterior probability distribu- 
tion for one or two parameters alone. This marginalized pos- 
terior is given by integrating over the other parameters. For 
example: 

P(9 u 9 2 \d) = J ni 3 d6iC0)P V n O M (3) 

where n is the number of parameters. We use the prior to 
incorporate non-CMB information suc h as that the r edshift 
of reionization must be greater than 6.0 (Becker et al. 2001). 

In the following subsections we discuss first how we evalu- 
ate the likelihood function at a single point and then how we 
evaluate it over a large parameter space and produce marginal- 
ized posterior distributions. 

3.1. Likelihood evaluation 

The first step to likelihood evaluation is to calculate the an- 
gular power spectrum, C/, for the cosmological model. To do 
this we use CMBfast. Despite its speed advantage s we do not 
use the Davis Anisotropy Shortcut (DASh; Kanlingh at et alJ 



CHU & KNOX 



3 



since the split loi, modifications were easier with 
CMBfast than with DASh. 

Once C\ is calculated we evaluate the likelihood given the 
WMAP data with the subroutine available at the LAMBDA 1 
data archive. For CBI and ACBAR we use the offset log - 
normal approximation of the likelihood ( B ond et ail 2000). 
The likelihood given al l these data to gether (referred to as 
the WMAPext dataset in Sper gel et all (120031) ') is given by the 
product of the individual likelihoods. 

We do not use th e most recent releas e of CBI data 
jReadhead et alJl2004 . nor the VSA d ata (IDickinson et all 
120041) . We have used the older release (Pearson et al. 2003) 
for ease of compari son with results in I Snergel et al. (2003). 
The new CBI data (Read head et alJl2004l) and the VSA data 
are consistent with the old CBI data, WM AP and ACBAR 
( IDickinson et al 120041: iRebolo et al.l2004l) . 

3.2. Exploring the Parameter Space 

We explore several parameter spaces. What we refer to as 
the six-parameter model has six cosmological parameters (the 
baryon density, the cold dark matter density, the scalar pri- 
mordial power spectrum amplitude at k = 0.05 Mpc -1 , n$, the 
redshift of reionization and S1a) and a calibration parameter 
for each of CBI and ACBAR. The split model is the same 
except the baryon density, lu/,, is replaced with uji, e and u>b p . 

We explore the parameter space by producing a Monte 
Carlo Markov Chain (M CMC) via the Me t ropoli s-Hastings 
algorithm as described in Christensen et al. (2001). Our pro- 
cedure is the same as in IChu et alJ II2003) except for some 
changes to the adaptive phase of the sampling, during which 
the generating function is determined. Our WMAPext chain 
for the six-parameter model space has 80,000 samples and the 
WMAPext chain for the seven-parameter split model space 
has 130,000 samples. 

4. RESULTS 

We now examine how the likelihood function changes as 
the parameter space is expanded to non-zero u>b p —u)be- 

4.1. Constraints on ujbe and Wbp 

In Fig. 2we see the resulting constraints on u>b (assuming 
no splitting), Lo be and iOb P - Note that, without splitting, we 
find ujb = 0.0229 ± .0012 which reproduces the result from 
Sperg el et al.1 (|2003) of 0.023 ± 0.001 for the same model 
space and dataset. 

For the unsplit parameter there is a tension between increas- 
ing the damping length (which occurs by lowering uib) and 
increasing inertial effects (which occurs by raising uob). This 
tension only becomes evident once we split the parameter and 
see that Lo be drifts downward and Lo bp drifts slightly upward. 

In Fig. [3] we see why the data tend to favor lower to be- it 
improves the fit in the damping tail region. Model fits with 
the six-parameter model lead to excess model power at higher 
/. Allowing ujbe to vary independently from u>b p gives the nec- 
essary freedom to eliminate this excess and improve the over 
all fit. The tendency for u>b p to drift upward, seen in Fig. 2, 
has indirect causes as we explain in the next subsection. 

The excess small-scale power of the six-parameter model 
fits is also what causes the slight p reference of the data for 
dns/dlnk < 0. iSpergel et all (120031) find for WMAPext that 

1 The Legacy Archive for Microwave Background Data Analysis can be 
found at http://lambda.gsfc.nasa.gov/ 



dns/dlnk = -0.055±0.038. This extension by splitting ui/, ac- 
tually leads to somewhat better agreement with the data than 
the extension to dns/dlnk 0. We find the maximum likeli- 
hood improves by a factor of 8.2 with the extension to split ujb 
and 2.8 with the extension to dns/dlnk^O. 

In the right panel of Fig. 2 we see the constraints in the 
0Jb P -i>Jbe plane. The straight line is the physical subspace, 
u>bp = Ube- How significant is the deviation from the physi- 
cal subspace? We find u) bp -w be = 0.0088 ± 0.0039, a 2.3<r 
difference from zero. The difference from zero could simply 
be a statistical fluctuation, it could be caused by systematic 
error in one or more of the experiments, or it could be that the 
six-parameter model does not adequately describe reality. The 
discrepancy is not strong enough to rule out the first option. 

4.2. Constraints on other parameters 

Splitting LOb alters the constraints on other parameters too. 
To gain a better understanding of the model dependence we 
now examine how constraints on r, n$ and uj„, change. 

In the six-parameter model there is a tension between the 
value of t that best fits the Cj T data (r = 0) and the larger 
value which best fits the Cj E data. Allowing extra freedom in 
the damping tail, by splitting ujb, makes it possible for a higher 
r to be consistent with the Cj T data as can be seen in Fig. 4. 
Increasing n$ can compensate for increased r, but too large an 
n$ leads to too much power on small scales. Extra damping 
from decreased ujb e compensates for this excess power. Fur- 
ther, increased n$ increases the ratio of the 2nd peak height to 
first peak height, an d ther efore ujb p increases to compensate. 

We also see in Fig. l4.2l that the likelihood fu nction of u> m has 
broadened some toward higher values. Spergel et al. (2003) 
claim that the constraint on u> m is largely coming from the 
rise to the first peak and how the shape of the spectrum in 
this region is influenced by the early integrated Sachs-Wolfe 
(ISW) effect. The amount of fluctuation power from the early 
ISW effect depends on the ratio of matter energy density to 
radiation energy density at last scattering. We see from our 
split model chains that models with higher z* also tend to have 
higher oj m , as they must to keep p m j p z at z* nearly fixed. 

4.3. Benefits of measuring the damping tail 

Measuring the damping tail region of the spectrum to high 
precision will be very valuable. In the context of the six pa- 
rameter model it will allow for much tighter constraints on 
LOb, oj,„ and ns- Whereas WMAP will measur e out to about 
/ = 500 with cosmic variance precision, Planck ( Taube rl2001l) 
will measure out to about I = 2000 with cosmic variance pre- 
cision, reducing the errors on the above paramete rs by factors 
of 6, 5 and 6 respectively (Eis enstein et alJl9 99). 

Perhaps more importantly, measuring the damping tail will 
allow more stringent tests of the model. The six-parameter 
model, calibrated with measurements at larger angular scales, 
makes tight predictions for the damping tail region. For ex- 
ample, the damping scale Id, is very tightly constrained with 
the WMAP ext d ata given the six-parameter model, as can be 
seen in Fig. 14. 21 But these tight constraints are not because we 
have measured the damping tail well. They are due to the fact 
that the parameters controlling the damping tail region can be 
determined well at larger angular scales. 

With the split model we weaken this connection between 
the acoustic region and the damping region. As a resul t, th e 
constraints on Id weaken considerably as seen in Fig. 14.21 . 
Future high precision measurements of the damping tail will 
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FIG. 2. — Left panel: One-dimensional marginalized likelihoods of u>t (assuming no parameter splitting), W] x and u>tp. Right panel: Two-dimensional 
marginalized likelihood of uji, p vs. tote- Dark (light) shaded region is where the likelihood is down by less than a factor of exp(-2.3/2) (exp(— 6.17/2)) from its 
maximum value. The line is uj^p = u>be- 



1= 



1000 



■a 500 



CO 

'm 
K 




-500 



1000 



1500 



1 



FIG. 3. — Top panel: CMB constraints on the power spectrum, the fiducial 
six-parameter model, and an improved fit that comes from allowing uj[, p ^ 
uite ■ Bottom panel: Same as above, but with the fiducial six-parameter model 
subtracted. 



be able to determine lo with high accuracy. Then the predic- 
tion for la, given the six-parameter model and acoustic region 
data, can be compared with Id inferred from the damping re- 



gion, allowing for a strong test of the model. 

Current measurements of the damping tail are already play- 
ing an important role in testing the six-parameter model. 
ACBAR and CBI do reduce the allowed region of the six- 
dimensional parameter space some, but mostly serve to 
confirm the WMAP predictions (albeit with slightly lower 
power). Once we allow for the parameter split, ACBAR 
and CBI place significant extra constraints on the parame- 
ter space. For example, constraints on the split u>i, model 
from the WMAP data mean that Ci 40 o = 805 ± 169 (where 
C t = /(/+ 1)C;/(2tt)). Inclusion of the ACBAR and CBI data 
change this to 695 ± 79, reducing the uncertainty by a factor 
of two. In contrast, without the ojj, split, adding ACBAR and 
CBI only reduces the C1400 error by 30%. 

5. CONCLUSIONS 

We have introduced the creation of cosmological meta- 
parameters as a tool for exploring the powerful, compli- 
cated, model-dependent constraints on cosmological parame- 
ters possible with measurements of CMB anisotropy. Such an 
exploration is useful for gaining a physical understanding of 
the origin of the constraints and for testing the consistency of 
the model. We see, as expected, uif, is mostly constrained via 
the observable consequences of its effect on the inertia of the 
pre-recombination plasma. Electron-scattering effects play a 
sub-dominant role. Determinations of cjj from the two differ- 
ent effects differ by 2.3a; i.e., they are marginally consistent 
with each other. 

It is amazing that a model with only six parameters can 
provide such a good fit to the WMAPext dataset. With the as- 
sumption of this model, fairly tight constraints are possible on 
param eters such as n* and u)/,. We see though, as others have 
seen llTe ffmark et al]l2"o"04l . that extending the model space 
by just one parameter can greatly broaden the constraints on 
other parameters. Interpretation of CMB constraints on cos- 
mological parameters must be done with care. 
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Fig. 4. — Two- and one-dimensional marginalized likelihoods. Left panel: The dashed (solid) line is where the likelihood of n s and r given the WMAPext 
dataset assuming the six-parameter model is down by a factor of exp(— 2.3/2) (exp(— 6.17/2)) from its maximum value. Shading indicates the same regions, but 
for the split model space. Middle and right panels: One-dimensional marginalized likelihoods of u) m and Id given the WMAPext dataset. Solid lines are for the 
six-parameter model space and dashed lines are for the split model space. 



In particular we have seen that the modeling of processes 
that affect the damping tail region of the spectrum can greatly 
loosen bounds on n$ and r. Other (physical) parameter exten- 
sions that would affect the damping tail reg ion are ru nning, 
and a time-varying fine structure constant (Kaplingh at et al.l 
1999). Measurement of the damping tail with great precision, 
such as will be done by Planck, will dramatically decrease 
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